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Abstract 

The scaling properties of the time series of asset prices and trading volumes of stock 
markets are analysed. It is shown that similarly to the asset prices, the trading 
volume data obey multi-scaling length-distribution of low-variability periods. In 
the case of asset prices, such scaling behaviour can be used for risk forecasts: the 
probability of observing next day a large price movement is (super-universally) 
inversely proportional to the length of the ongoing low-variability period. Finally, a 
method is devised for a multi- factor scaling analysis. We apply the simplest, two- 
factor model to equity index and trading volume time series. 
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1 Introduction 

Predicting future developments of financial asset prices is an extremely difficult 
task. Certain problems, such as predicting the direction of movements, are 
nearly unsolvable. Meanwhile, certain aspects can be easily predicted. For 
instance, one can be sure that the prices will always fluctuate intermittently, 
largely due to people who believe they have found the winning algorithm, the 
philosopher's stone. 

The analysis of the historic charts of price dynamics is referred to as techni- 
cal analysis [1]. The essence of the technical analysis is the hypothesis that 
the patterns of historic data can forecast the future price movements. On the 
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other hand, the efficient market hypothesis (c.f. [2]) states that security prices 
reflect fully all available information. In its century-long history dated back 
to the work of Bachelier [3], the financial analysis has made use of the both 
approaches. The random walk hypothesis, standard deviation and the corre- 
lations between securities' returns are the cornerstones of seminal papers in 
financial analysis [4,5,6]. Econophysics has introduced various more elaborate, 
mostly non-linear tools of analysis (c.f. [7,8] and references therein). Here we 
extend a recently developed technique of length-distribution of low- variability 
periods [9]. 

Strongly non-linear systems, such as turbulent fluids and plasmas, granular 
media, biological and economical systems, etc., are typically characterized by 
scale-invariance and scaling laws. So, it should be not surprising that such 
seemingly different disciplines like turbulence studies, biological physics, and 
econophysics have many common tools of data analysis. For example, power- 
law distributions were first used in economics [10] in the end of 19th century, 
and later found in a wide variety of systems (often referred to as the Zipf's 
law), c.f. [11,12,13,14]. Similarly, diverse systems are known to generate sig- 
nals, which are self-affine, and hence, can be characterized by the Hurst expo- 
nent, c.f. [15]. Further, the stable Levy distributions (c.f. [16]) have been found 
to be relevant to all the mention disciplines; about truncated Levy flights in 
econophysics, c.f. [17]. Finally, multi-fractal formalism is by far the most pop- 
ular tool for scale-invariant analysis of intermittent time-series, c.f. [18,19]. It 
should be noted, however, that in many cases (including econophysical appli- 
cations), there is no profound understanding of the origin of multi-fractality, 
and hence, the multi-fractal analysis is not necessarily the optimal method. 
Indeed, multi-fractal formalism has been devised in the context of turbulence, 
and is specifically suited for systems with random multiplicative cascades [20]. 
Meanwhile, in the case of such time-series as stock prices or heart rate vari- 
ability signal, the presence of multiplicative cascades is not evident. 

Therefore, there is a clear need for a deeper understanding of the charac- 
ter of intermittency in the case of financial time-series. This problem can 
be approached by studying new independent and/or more general methods 
of scale-invariant analysis. For instance, in the case of heart rate variability, 
it has been found that the whole time-period can be clusterised into self- 
similarly distributed segments of approximately constant heart rate [21]. In 
order to address this problem, a new method has been suggested recently [9] 
which is based on the analysis of the distribution of low-variability periods. 
The low-variability period is defined as a time period with maximal length 
where consecutive relative changes in realizations of the time series are less 
than given threshold 5. Note that in addition to the financial assets, the low- 
variability period analysis has been applied to biological systems [22,23]. It 
has been shown [9] that in the case of the multi-affine time series, the cu- 
mulative distribution function of the low-variability periods is in the form of 
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a (multi-scaling) power-law. Since the opposite is not necessarily true, the 
low-variability period analysis is, indeed, a more universal method than the 
multi-affine analysis. 

Even if the time-series is actually multi-affine, the low-variability period anal- 
ysis can be still useful, because (a) the power-law exponent is related to the 
multi-fractal dimension; this circumstance provides an easy method for check- 
ing the assumption of multi-affinity; (b ) low- variability periods provide higher 
time-resolution of time series analysis [9] . 

This paper serves three main purposes. First, we are going to apply the method 
of the analysis of low- variability periods to the data of trading volumes. This 
analysis is motivated as follows. Similarly to the asset prices, the trading vol- 
umes are known to fluctuate intermittently, c.f . [24,25] . According to the Man- 
delbrot's model of stock prices as a fractional Brownian motion in multi-affine 
trading time [18], one could expect that the time series of trading volumes are 
multi-affine. The analysis of the length-distribution of low-variability periods 
of trading volumes serves as a test of this model. Second, we discuss the con- 
sequences of the presence of power laws of low-variability periods. Third, an 
attempt is made to generalize the method to multivariate timc-scrics (e.g. a 
stock price together with the trading volume). It should be noted that in the 
case of multi-affine analysis, there is no simply interpretable way for such a 
generalization. 



2 Scaling of trading volumes 

We start with a brief description of the method devised in Ref. [9,22], A low- 
variability period is defined as such a contiguous time-interval Ti — [ti, U -\- Zj], 
which satisfies the conditions 



and \ V{t) I (yif))^ — 1| > 5 for t = — 1, -|- 1; it is assumed that the data 
sampling interval serves as the time unit. The average of the trading volumes 
V{i) is taken over a time window of length r: 



Therefore, wc have two control parameters: (i) r — the length of the averaging 
window, and (ii) b — the variability threshold. Further, the cumulative length- 
distribution function of low- variability periods R{n) is introduced: R{n) is the 
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number of the low- variability periods of length k > n. We speak about multi- 
scaling behaviour, if the following power-law is observed: 



where a{5,T) is a scaling exponent and Rq is a constant. 

It has been shown [9] that in the case of multi-affine time-series, the low- 
variability periods do follow a multi-scaling distribution law: 



where f{h) denotes the Holder multifractal spectrum of the local Hurst expo- 
nents h. Meanwhile, the presence of a scaling law (3) does not necessarily imply 
multifractality. Indeed, note that the Holder exponent f{h) cannot exceed the 
topological dimensionality (one), therefore the values a > 1 are not related 
to multifractality. It should be stressed that unequality r) > 1 does nei- 
ther imply the lack of multifractality: in the case of multi-fractal timc-scrics, 
a{6,r) > 1 can (in principle) be observed for log^5 < ho, where ho is the 
dominant local Hurst exponent, f{ho) — 1. For this range of parameters, the 
unequality a{S, r) > 1 must be satisfied (assuming that « is a monotonous 
function of d). Finally, if there is no data collapse a{6,T) = f{log^6) for 
o^i^, t) < 1, the underlying time-series is certainly non- multifractal. In the 
case of currency rates, a reasonable data collapse [according to Eq. (4)] has 
been observed. In this Section, the daily trading volume data of various stock 
indices are tested for multi-scaling and multi-fractality. 

2.1 The data and the calculations 

The stock market volumes are measured in the amount of shares traded on 
the exchange. The data used in the analysis represent the daily closing prices 
and trading volumes and it is described in detail in the Table 1. 

If the cumulative distribution function R{n) is plotted against n in log-log 
scale, the power-law Eq. (3) corresponds to a straight line. The scaling expo- 
nent a{d,T) was found as the slope of the linear trend line in log-log space, 
using the least-square fit method. 

The error of a{d, t) was estimated as follows: The least-squares fitted trend-line 
was found as described above, except that the slope a was not optimised, i.e. it 
was considered as a fixed parameter. Further, the sum of squared residuals r{a) 
was calculated as a function of a. The error estimate was found as e = (a' — a). 




(3) 



a{,5,T) = /(log^5). 
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Table 1 



The data used in volume analysis 



Abbr 


Description 


Calendar Period 


# of data 


SPX 


Standard & Poor's 500 Index 


04/01/93 - 13/09/04 


2947 


DAX 


The German Stock Index 


04/01/93 - 13/09/04 


2950 


NKY 


Nikkei 225 Stock Average 


04/01/93 - 14/09/04 


2879 


MXEA 


The MSCI Europe, Australasia 
and Far East Index 


04/06/01 - 13/09/04 


832 


CAC 


CAC-40 Index of Paris Bourse 


04/01/93 - 13/09/04 


2953 


UKX 


FTSE 100 Index 


04/01/93 - 13/09/04 


2953 


MXWO 


MSCI World Index 


04/01/01 - 13/09/04 


809 


INDU 


Dow Jones Industrial Average 


04/01/93 - 13/09/04 


2941 


TALSE 


Tallinn Stock Exchange Index 


25/02/03 - 13/09/04 


392 


RTS 


Russian Trading System Index 


03/07/01 - 13/09/04 


794 


WIG 


Warsaw General Index 


23/05/01 - 13/09/04 


836 


BUX 


Budapest Stock Exchange Index 


15/09/97 - 13/09/04 


1747 



where a is the least-squares fitted value of the slope, and a' satisfies the 
condition r{a') = 2r{a). 



We used S = [5%, 10%, 70%] and r = [2, 3, 5, 10] days as input parameters 
for all of the time series described in Table 1. Total 650 calculations were 
carried out and in most cases the least-square fit provided good results. The 
quality of the fit is measured by the coefficients that are plotted into the 
histogram in Fig 1. The scaling exponents (found with r = 10 days) with error 
estimates are presented in Table 2. 

2.2 Scaling properties of volume data and discussions 

As seen from Fig. 1 and Table 2, the low- variability periods of the daily vol- 
ume time series follow reasonably well a multi-scaling behaviour (similarly to 
what has been observed in the case of stock prices [9]). In Fig. 2, the scaling 
exponent a{S,T) is plotted against the logrS for DAX and SPX time series 
with T = 3, 5, 10 days. We can see that there is a departure from the multi- 
affine expectation — there is no data collapse at range o<l [which would 
have been corresponding to Eq. (4)]. This observation is interpreted as follows. 
Here, each data point describes the amount of shares traded in respective stock 
exchange in a calendar day. The condition a < 1 is satisfied only for very high 
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Table 2 



The values of the scaling exponent a{6, t = lOdays) in the case of daily trading 
volumes for various stock indices and different values of the parameter 6 
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Fig. 1. Histogram of R-squared coefficients based on the regression analysis of de- 
termining scafing exponent a((5, r) 
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thresholds {5 > 50%). Typically, there are very few periods, when trading 
volumes would fluctuate with such a large amplitude. Therefore, the number 
of valid data points is very low. As a consequence, the scaling range is narrow, 
and the formally calculated scaling exponent is non-stationary (depends on 
the particular realisation of the time-series). This behaviour is different from 
what has been observed for currency time series (when the results a <1 were 
more stationary, with a reasonable data collapse in a — log^ (5-plot). 

So, we can conclude that while the trading volume data are intermittent and 
fluctuations are scale-invariant [described by Eq. (3)], the degree of intermit- 
tency is lower than in the case of currency fluctuations: the scaling exponent 
values tend to be always larger than one. Respectively, the multifractal pattern 
[like described by Eq. (4)] is not observable. 

2.3 Consequences of the scaling behaviour 

In the case of financial time-scries, a crucial question is how to make forecasts. 
The most useful kind of forecasts would give predictions of the direction of 
future prices. However, these forecasts will be always very unreliable, the ef- 
ficient market hypothesis denies such a possibility entirely (c.f. second law of 
thermodynamics). Meanwhile, the risk-related forecasts are completely possi- 
ble (and may prove to be useful). Therefore, a natural question arises: what 
are the risk-prediction-related consequences of the presence of the power-law 
distribution of low- variability periods? More specifically, suppose we have had 
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Fig. 2. Scaling exponent a(5, r) plotted against loQrS for SPX and DAX time series 
for T = 3, 5, 10 days 
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a low- variability period, which had lasted n days (and still goes on). Can we 
say something about the possibility of today being the last day of this "silent" 
period? In other words, what is the probability p{n) that the tomorrow's 
movement exceeds our pre-fixed threshold 5? 

Apparently, this probability is given by the ratio of (a) the number of those 
low- variability periods, the length of which is exactly n, Na = R{n)—R{n + 1), 
and (b) the number of those low- variability periods, which have length m > n, 
Nf, = R{n). If n is large, the difference R{n) — R{n + 1) can be calculated 
approximately as — Upon applying Eq. (3) we arrive at Na ~ ai?o^~"~^ 
and Rqu'"'. Bearing in mind that p{n) 

written as 



Na/Nb, the final result is 



an 



(5) 



Therefore, we have shown that the very presence of a power law of the low- 
variability periods has an interesting consequence: the probability that the 
tomorrow's price movement will be larger than the movements of n pre- 
ceding days, is inversely proportional to n. The predicted power law expo- 
nent is independent of the scaling exponent r), i.e. we are dealing with 
super-universality. This super-universality appears to be related to the super- 
universality of the scaling of direct avalanches in self-organised critical systems 
[26] . It allows us to make every day a series of forecasts: by scanning the values 
of the threshold parameter 6, we find the length ns of the current (still ongoing) 
low-variability period. Then, the probability that the tomorrow's movement 
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p{S) remains below the threshold 6 is estimated as p{5) ~ a{6, r = 2 days)n^''". 
Note that the prefactor a{6, r = 2 days) ~ 1 can be dropped, but will improve 
the forecast; r = 2 days means that only the prices of two subsequent days 
are compared. 



3 Two-factor model of scale invariance in financial time series 

In the technical analysis, various forms of multi-factor models have been used 
for some time [1]. In particular, it has been conjectured that in the case of trad- 
ing volumes and price fluctuations, the higher-than average trading volumes 
generally "confirm" the price trend. However, according to our best knowl- 
edge, scale-invariant methods have not yet been applied to the multi-factor 
studies. 

The low-variability period analysis provides simply interpretable way of gen- 
eralisation of multi-signal time series analysis. With the conventional multi- 
affine analysis, the interpretation of the results is significantly harder. Indeed, 
consider the wavelet transform w{a,b), where a denotes a coordinate and b 
stands for the wavelet width. For a multi-affine signal, the following scaling 
law is expected: {w{a,by) oc 6^''. For a bi-variate analysis of two signals, this 
equation can be generalized as {wi{a.bY ■ W2{a,bY) oc b'^^-p. Here, Wi and W2 
denote the wavelet transform of the two signals. For two dependent signals, 
{wi{a, by • ^2(0, by) 7^ {wi{a, by) ■ (^2(0, 6)^); however, there is no clear way 
of interpreting the features of the scaling exponent Tg^p. Therefore, in this sec- 
tion, the analysis of low- variability periods is generalised and the multi-factor 
model is proposed. In particular, we address the above mentioned statement 
about the relationship between trading volumes and stock prices. 

In our previous calculations, the low-variability periods are defined either via 
the condition (1), or via the same unequality, with trading volume V{t) being 
substituted by stock price (or index value) p{t) [9]: 

\p{t)/{p{t))^-l\<Sp for teT,. (6) 

Hereinafter we refer to the application of Eq. (6) as the usage of the single- 
factor price model or Method (for brevity). Similarly, the condition (1) cor- 
responds to the single-factor volume model or Method (Volume). There are 
two ways to generalize the concept of low- variability into multi-factor model: 

(1) Low-variability persists if both condition (1) and (6) are satisfied; 

(2) Low-variability persists if any of the conditions (1) and (6) are satisfied. 

It is clear that the number of low-variability periods is larger by using the 
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latter option. In fact, using the first option led us often to a very small number 
of low-variability periods (this, of course, is related to the limited length of 
the time-series). Therefore, here we present only the results corresponding to 
the second definition. Note that for multi-factor models with more than two 
inputs, the conditions similar to Eq. (6) can be combined with any set of 
logical operators "and" and "or" . 

In what follows we use several definitions of the two-factor low-variability 
periods (the apphcation of these definitions will be referred to as Method 1 - 
Method 3, respectively): 

- 1| < ^) Vd^W/ (^W)r - 1| < ^v) for t e T,. (7) 

m)/ m)r - ii < ^v) \J(y{i)i (^w)r - 1 < for t e (s) 
iw)/ m)r - ii < ^) yim/ (^w). - 1 > -^.) for t e r,. (9) 

3.1 Asymptotic behaviour 

Equation (7) is symmetric with respect to the price and volume conditions. 
The low-variability period is terminated as soon as the relative price change 
exceeds (by modulus) Sp, and the relative volume change exceeds (by modulus) 

If the condition (7) is applied, the following asymptotic laws are expected: 

a{Sp,6v,T) — >• a{5p,T) when 6v — >■ 0, 

and 

a{5p, 5v, t) — > a{5v, t) when Sp — > 0. 

Such a behaviour is, indeed, observed for the time series of DAX and SPX. 
The data used in the calculations is described in Table 1. In Figures 3 and 4 
the scaling exponents a{dp, 6v, t) are displayed for DAX and SPX time series 
by using r = 10 days. iFiom Fig. 3, it is seen that for low values of 6v, the 
low-variability condition for volume is almost always violated. Therefore, the 
scaling exponent is defined almost solely with the price parameter (the curves 
with Sy — and Sy = 0.5 are very close to each other). On the other hand, for 
high values of Sy, the volume condition is typically satisfied, and the horizontal 
lines (in a-6p space) denote a low dependence on the price parameter. These 
arguments apply to the dependence of the exponent a{dp,Sv,T) on the price 
threshold parameters, as well, see Fig. 4. 
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Fig. 3. The convergence of Method 1 Method 0. The scahng exponents a{5p, 6v, t) 
are calculated for DAX time series by using conditions (6) and (7) with r = 10 days 
and 5y = [50%, 25%, 15%, 5%, 0%] 
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Fig. 4. The convergence of Method 1 — > Method (Volume). The scaling exponents 
a((5p,5y,r) are calculated for SPX time series by using conditions (1) and (7) with 
r = 10 days and 5p = [2.5%, 1.9%, 1.0%, 0.1%, 0.0%] 




5.0 



25.0 



45.0 



65.0 



5v(%) 



MO(Volume) 


-^Ml 


5p=0.0% 


^Ml 


5p=0.1% 


-B-Ml 


dp=1.0% 


^|«-M1 


d>1.9% 


^<-Ml 


(5p=2.5% 



3.2 Methods 2 and 3: differences between volume spike and volume squeeze 



According to Eq. (8), the low- variability period is terminated when price 
change (rise or drop) is significant (i.e. larger than the threshold parame- 
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ter), and the volume increases faster than the threshold. Equation (9) repre- 
sents a definition, opposite to Eq. (8): the low- variability period is terminated 
when the price change exceeds the threshold parameter, and the volume de- 
creases faster than the threshold. These definitions are useful for studying the 
asymmetry between the volume rise and drop: if the multi-scaling exponent 
a{5p, 5v, t) turns out to be different for Methods 2 and 3, there must be an 
asymmetry between those volume spike and squeeze events, which are accom- 
panied by a large price variability. Indeed, the price condition in Eqns (8) and 
(9) is the same; so, the differences in the scaling exponent a{Sp,Sv,T) must 
be due to the different effect of the volume condition (8) and (9). 

Let us refer back to Method 1 [Eq. (7)]. With this method, the events termi- 
nating the low-variability periods represent a superposition of the respective 
events for the Methods 2 and 3 [Eqns (8) and (9)]. So, if the scaling exponents 
calculated according to Method 2 are very similar to the ones calculated ac- 
cording to Method 1, we can conclude, that the amount of the low- variability 
periods defined by Method 3 is insignificant. The same holds true if the scaling 
exponents of Method 3 tend to be similar to the ones of Method 1. 

The multi-factor scaling exponents a{Sp,6v,r) are found for the DAX data 
with r = 10 days. In Fig. 5, the scaling exponents a{Sp,Sv,T = 10) of DAX 
index are plotted against logr=ioSp with (a) 6v = 5%, (h) 5y = 15%, (c) 
Sv = 25% and (d) Sy = 50% . The Methods 0-3 [Eqns (6)-(9)] are used for 
definition of low-variability periods. 



3.3 Discussion of the results 

By giving a small value to the thresholds 6v and 5p in any of Eqns (6)-(9), the 
number of longer low- variability periods becomes small. Therefore, we have a 
large number of short low- variability periods (compared to high S values); this 
leads to high values of the scaling exponent. In Fig. 5a, the threshold value of 
dv is set to very low level of 5%. From Eq. (7) it can be seen that then the 
two-factor model leads to results which are very similar to the single-factor 
(i.e. price factor) model. Likewise, the larger the volume threshold Sy, the 
larger the difference between the a-5p-curves of the single- and multi-factor 
models (at Sy — 50%, the curves are rather dissimilar). 

An important issue is the difference between the curves of Methods 2 and 3. 
One can notice that there is a different behaviour at the small values of the 

parameter Sp, an evidence of the asymmetry between volume rise and drop. 
One can also notice that the scaling exponents calculated according to the 
Method 3 [Eq. (9)] are lower than the ones calculated according to Methods 
1 and 2 [Eqns (7) and (8)]. Meanwhile, the difference between the outcomes 
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Fig. 5. Multi-affine scaling of the two- factor model. The scaling exponents 
a{6p,6v,T) of DAX time series are plotted against logr=io6p with (a) 5y = 5%, 
(b)6v = 15%, (c) 6v = 25% and (d) 6v = 50% 
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of the Methods 1 and 2 [Eqns (7) and (8)] is minor. Therefore, we conclude, 
that high price variability is typically accompanied by increasing volume. This 
conclusion is independent of the price/volume pre-history (i.e. is valid both 
for short and long low- variability periods). 

In this paper, we have not analysed the problem of higher-rank multi-factor 
models. However, this can be useful for e.g. multi-stock data analysis, where 
each stock price provides an independent input stream. This situation will be 
addressed in further studies. 



4 Conclusion 



The concept of low- variability periods has been proven to be useful for various 
econophysical issues (not just limited to the scope of stock prices/indices and 
currency exchange rates). So, we found that the time series of stock trading 
volumes obey multi-scaling properties, similarly like the price data. However, 
while the multi-scaling exponent of the price time series follows a pattern, 
characteristic to the multi-affine data, in the case of trading volumes, there is 
a clear departure from that pattern (one can say that the fluctuations are less 
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intermittent) . 

Further, we have shown that the presence of the multi-scahng distribution 
of the length-distribution of the low-variabihty periods gives rise to a super- 
universal scaling law for the probability of observing next day a large price 
movement. This probability is inversely proportional to the length of the on- 
going low-variability period, a fact which can be used for risk forecasts. 

Finally, the multi-factor model is proposed for time series analysis. In this pa- 
per, only the simplest two-factor model is described and applied to stock price 
and volume data. The low-variability periods of multi-variate time series can 
be defined in different ways; for instance, the threshold conditions apphed to 
the single data streams can be combined by logical "and" , as well as by logical 
"or" . In the our case of price and volume data, three different definitions of 
low- variability periods have been applied (in order to study the asymmetry be- 
tween the volume rise and drop, we have also applied sign-dependant threshold 
conditions). This analysis led us to the conclusion that high price variability is 
typically accompanied by increasing trade volume, independently of the prior 
events of the market. In the fight of this observation, the common thesis of 
technical analysis, "increased trading volumes confirms the price trend" , be- 
comes less useful. Indeed, most of the significant price jumps are accompanied 
by increased trading volumes; so, almost all the "price trends" pretend to be 
"confirmed" . 
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